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Solitons are of fundamental importance in photonics due to applications in optical data transmission 
and also as a tool for investigating novel phenomena ranging from light generation at new frequencies 
and wave-trapping to rogue waves. Solitons are also relativistic scatterers: they generate refractive- 
index perturbations moving at the speed of light. Here we found that such perturbations scatter 
light in an unusual way: they amplify light by the mixing of positive and negative frequencies, as we 
describe using a first Born approximation and numerical simulations. The simplest scenario in which 
these effects may be observed is within the initial stages of optical soliton propagation: a steep shock 
front develops that may efficiently scatter a second, weaker probe pulse into relatively intense positive 
and negative frequency modes with amplification at the expense of the soliton. Our results show a 
novel all-optical amplification scheme that relies on relativistic scattering. 

If compared to the well-developed field of traditional light scattering in which the medium is at rest, little attention 
has been devoted to the physics of scattering from a moving medium, in particular from a relativistically moving 
medium. Here we consider the remarkable ability of solitons to generate a co-propagating refractive index inhomo- 
geneity that propagates at relativistic speeds. The basics of scattering from a time-changing boundary were discussed 
in detail by Mendonga and co-workers (see e.g. [1] and references therein). Examples of such "time refraction" have 
been predicted and observed from a moving plasma front [1-3] and in waveguide structures [4-6]. Recently, the 
nonlinear Kerr effect, i.e. the local increase of the medium refractive index induced by an intense laser pulse [7], was 
proposed to induce a moving refractive index inhomogeneity within a dispersive medium such as an optical fibre [8] . 
The laser pulse induced relativistic inhomogeneity (RI) was then described in terms of a flowing medium in which the 
analogue of an event horizon may form [8] and applications such a optical transistors have been proposed [9]. 
Intense laser pulses are also known to scatter from the self-induced travelling RI: this self-scattering process leads 
to the resonant transfer of energy from the laser pulse to a significantly blue-shifted peak, often referred to as reso- 
nant radiation (RR) or "optical Cherenkov" radiation [10-14]. A recent discovery highlighted an additional scattered 
mode, further blue-shifted with respect to the RR, identified as a mode excited on the negative frequency branch of 
the medium dispersion relation and therefore named "negative resonant radiation" (NRR) [2]. 

By applying a first Born approximation analysis, supported by numerical simulations, we show that the scattering of 
light by a soliton induced RI reveals a novel all-optical amplification mechanism. A schematic representation of the 
specific process we are considering is shown in Fig. 1: an incoming laser pulse, IN, interacts with a co-propagating 
inhomogeneity, or scatterer. The momentum conservation law that governs the scattering process predicts that light 
may resonantly scatter into two output modes, RR and NRR (see Additional Material). In the laboratory reference 
frame, both of these modes will have positive frequencies while in the reference frame comoving with the scatterer, 
RR is positive and NRR is negative valued. The key point in the following is precisely the presence of these two 




FIG. 1: Schematic representation of scattering from a relativistic inhomogeneity: an incoming laser pulse, IN, interacts with a 
co-propagating RI, or scatterer, that transfers energy into two output modes, RR and NRR. 
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FIG. 2: Scattering from a "toy-model" RI in diamond: IN pulse at 2 ^m; super-Gaussian moving RI with m = 40 (3 fs rise- 
time), (Sno = 0.08. Envelope (a) and spectral (b) evolution along propagation, (c) Photon number evolution along propagation: 
the total sum (photon number amplification) increases with increasing RI steepness (increasing parameter m). (d) Comoving 
dispersion curve with indicated the allowed optical modes. 



modes: mixing between the RR and NRR modes leads to a novel amplification process. 



RESULTS 



First Born approximation. In scattering theory, the first Born approximation allows us to estimate the elements 
of the scattering matrix S, which connects asymptotic input and output states. Kolesik et al. have shown that the 
positive RR mode is captured in great detail by the Born scattering approach [3] . Here we extend this theory in order 
to also include the possibility of scattering to the NRR mode and find that the amplitude of each scattering channel 
is proportional to the amplitude of the Fourier component of the RI, R{uj) (see Additional Material): 

V - 

S{z,UJrr) W -i?(wRR - tJiN) (1) 

and, accounting also for the negative frequency branch, we similarly obtain: 

V " 

'S'(2:,t^NRR) ~ -^R{^NRR + ^m), (2) 

where v is the velocity of the travelling RI. These relations (1) and (2) state that energy may be transferred to two 
modes, RR and NRR, and in order to do so the spectrum of the scattering potential must have non-zero Fourier 
components at frequencies equal to the relative distances between cjin and the resonant frequencies. In other words, 
the RI must present a sufficiently steep gradient so as to effectively excite the desired output modes. We now observe 
that we may derive a photon number balance equation by generalizing the Manley-Rowe relation, adopted e.g. in 
nonlinear optics [7], to the case of a moving scatterer [4, 5]. We find that: 

|RRp - |NRR|2 = 1, (3) 

where |RRp and |NRRp are the photon numbers of the RR and NRR modes normalized to the input photon number, 
|INp. The negative sign in front of the |NRRp photon number is a direct consequence of the fact that the NRR-mode 
has negative frequency in the comoving reference frame (see Additional Material). So the difference between the 
normalized number of photons has to be equal to the photon number in the input mode. As a consequence, the 
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total output photon number, |RRp + |NRRp > 1, i.e. we have amphfication. The scattering process mediated 
by the traveUing RI will amplify photons as a result of the coupling between the positive and negative frequency modes. 

Numerical simulations. Wc verified our predictions by numerically simulating the one dimensional propagation 
of a scalar field propagating along the z direction [19] in a transparent and isotropic dielectric medium (diamond 
and fused silica) that is arranged to scatter from a co-propagating inhomogeneity (see Methods). We consider three 
separate cases: (i) the underlying physics arc first exemplified by taking a "toy-model" situation in which the RI is 
simply a variation in the linear refractive index; (ii) the RI is generated by an actual soliton-like laser pulse through 
the nonlinear Kerr effect and light from the soliton itself is self-scattered; (iii) the RI is generated by a soliton-like 
pulse and a second, very weak probe pulse is scattered by the soliton. 

(i) - ^'Toy-model" RI. The RI is simulated as a linear propagating refractive index inhomogeneity that moves 
with a speed v = 1.23 x 10^ m/s that is just slightly slower than the group velocity of an input probe pulse in the 
medium, Vg = 1.25 x 10* m/s. In this C3jSG, clS Si dielectric medium we choose diamond, which does not exhibit any 
resonances over a very broad bandwidth (from less than one Terahertz through to ultraviolet wavelengths) . The probe 
pulse is very weak so that no nonlinear effects are excited and the physics are dominated solely by linear scattering 
from the RI. Figures 2(a)-(b) show the evolution along propagation of the envelope (plotted in the frame moving at 
the RI velocity) and of the spectrum, respectively. We launch an input Gaussian probe pulse (carrier wavelength 2 
/um, pulse width 20 fs), together with a super-Gaussian shaped RI (see Methods), with amplitude 5no = 0.08 and 
super-Gaussian order m = 40 corresponding to a rising time ^ 3 fs; cr = 50 fs. This ensures a wide RI Fourier 
spectrum Sn{oj), so as to meet both conditions (1) and (2). After ~ 0.2 mm the probe pulse reaches the moving RI 
[whose position is highlighted in Fig. 2(a) by the vertical black dotted lines] and slows down due to the RI increase of 
refractive index. Scattering occurs towards two blue-shifted resonant modes, RR and NRR, which due to dispersion 
both travel slower than the RI. The medium dispersion relation is shown in Fig. 2(d) in {oj',u)) coordinates where 
the comoving frequency u)' = j{u — vk); k = u)n{u))/c is the medium dispersion relation and 7 = I/1/I — (w/c)^, 
being c the speed of light in vacuum. In the comoving frame, the absolute value of scattered mode frequencies arc 
determined by the input mode comoving frequency (horizontal solid and dashed lines, see Additional Material): the 
corresponding laboratory reference frame frequencies can be seen to be in perfect agreement with the numerically 
observed values. 

In order to verify the photon number amplification predicted by the Manley-Rowe relation (3), we evaluated the 
photon number evolution, as shown in Fig. 2(c) for two different RI gradients, obtained by varying the order of 
the super-Gaussian function that describes the RI, m = 40 and m = 200 (these correspond to rise times from the 
background refractive index to the maximum of the RI of 3 fs and 0.7 fs, respectively). Dotted, solid and dashed 
lines correspond to the normalized photon numbers in the IN mode, total output photons |RR|^ + |NRRp and output 
photon difference |RRp — |NRRp, respectively. As can be seen, the difference in photon numbers is conserved while 
the sum of photon numbers is larger than 1 (amplification) and increases with increasing RI steepness. This simplified 
"toy-model" therefore gives a direct confirmation of the predictions based on the first Born approximation model. 

(ii) - Soliton-induced RI. The RI is physically generated by an intense soliton-like laser pulse through the nonlinear 
Kerr effect, i.e. Sn = 712/, where / is the peak intensity and n2 the nonlinear Kerr coefficient (see Methods). The 
soliton may interact with the self-induced refractive index variation and scatter into new output modes. As a Kerr 
medium we choose common fused silica glass, whose dispersion relation is well known [20]. We launch a soliton-like 
mode in the anomalous dispersion region, with central wavelength 2 /xm, pulse width 7 fs and input intensity 10 
TW/cm^. Figures 7(a)-(b) show the envelope and spectral evolution, respectively, along propagation over a distance 
z = 0.4 cm. The input pulse initially self-steepens, forms a shock front and, at z ^ 0.3 mm, when conditions (l)-(2) 
are met, energy starts to resonantly transfer from the input mode towards two blue-shifted modes, RR and NRR, 
both travelling slower with respect to the input soliton. The soliton frequency is red-shifted due to a recoil effect [10] 
and subsequently emerges with a slightly lower propagation velocity. In Figure 7(c)-(d) we compare the nonlinear 
propagation with a numerical simulation in which the nonlinear source term has been replaced with a RI in the 
linear polarization with the same amplitude, steepness and propagation velocity of the soliton induced perturbation 
of Figs. 7(a)- (b). We observe a nearly perfect agreement with the nonlinear case confirming that soliton shedding of 
RR and NRR modes is nothing more than a specific realization of the more general relativistic scattering process. 
The "toy-model" RI scattering is indeed able to capture all the essential features of the resonant energy transfer, and 
coupling between positive and negative frequency modes. 

(iii) - Probe pulse scattering from a soliton-induced RI. Finally, we consider the case in which the soliton is accom- 
panied by a second, delayed probe pulse. This second pulse is much weaker and therefore does not form a soliton or 
excite any nonlinear Kerr effects. Figure 4 shows the results for a 1.9 /im central wavelength probe pulse, pulse width 
8 fs and input intensity 5 x 10^ W/cm^ (complete movie showing the full dynamics is also included). The probe pulse 
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FIG. 3: RI scattering dynamics in fused silica. Comparison between a soliton-induced (via the nonlinear Kerr effect) RI (a)-(b) 
and the "toy-model" RI with the same shape as the soliton-induced inhomogeneity, in a purely linear medium (c)-(d). Envelope 
(left hand side), and spectral (right hand side) evolution in propagation. 
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FIG. 4: Weak probe scattering in fused silica, from the soliton-induced RI of Fig. 7(a)-(b). Probe envelope (a) and spectral 
evolution in propagation (b). Normalised photon number evolution (c) for the probe (black solid curve) and for the soliton 
(blue dashed curve). A complete movie animation is included. 



is scattered by the soliton induced RI shown in Fig. 7(a)-(b). The initial probe pulse delay, — 15 fs, is adjusted 
such that it encounters the soliton when it forms the steepest shock-front and is then scattered relatively efficiently, 
simultaneously spectrally recoiling to a slightly red-shifted wavelength and emitting two blue-shifted RR and NRR 
modes. In Fig. 4(c) we show the overall photon number evolution for the probe pulse (black solid line) and for 
the soliton pulse (blue dashed curve). In order to highlight the cross-scattering dynamics, the two curves have been 
normalized at each propagation distance with respect to two independent simulations for the probe and soliton pulses 
alone. We thus note that the overall probe pulse photon number increases by nearly 1%, followed by some weaker 
oscillations that are originating from a remnant of the probe pulse that is trapped and thus continuously interacts with 
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the soliton [sec Fig. 4(a), trapped light is visible for z > 1.5 mm, around t ^ fs]. This clear photon number increase 
thus indicates that soliton-induced scattering in this regime, very differently from standard cross-phase modulation, 
may lead to true amplification of the weaker pulse. 

DISCUSSION 

Summarising, we have shown that a RI amplifies and scatters light to higher frequencies. Likewise, if the probe pulse 
were to be reduced to the level of quantum fluctuations, we may expect to see the RI excite the vacuum states. 
Bearing in mind that the motion of photons in the vicinity of a RI may also be described in terms of an effective 
curved space-time metric [21, 22], we expect that measurements of vacuum fluctuations excited by a RI would give 
an experimental direct window into the physics of the quantum vacuum in curved space-times [23-25] . 

These results therefore open the perspective for novel all-optical control schemes that may be implemented in a wide 
variety of geometries and applications together with a novel numerical and experimental approach for the investigation 
of fundamental phenomena. 

METHODS 

Numerical model. The code adopted for the numerical results presented in the main text is a ID Unidirectional 
Pulse Propagation Equation (UPPE) solver [19], where the spectral components Ei^{z) of the real field E{z,t), are 
obtained by Fourier transform and obey to the following model equation, in the pulse or inhomogeneity comoving 
frame: 

d,E^ - i[k,{w) - uj/v]E^ = . . P„, (4) 

2eoC^fcz(a;) 

being eo the vacuum permittivity, c the speed of light in vacuum, and P^^ the spectral component at frequency oj of 
the real-valued polarization as obtained by Foiiricr transform. 

Dispersion of the medium is described by kz{oj) and the shift —oj/v indicates that Eq. (4) is solved in the local frame 
of the scatterer. Below we report the "linear" (i) and "nonlinear" (ii) implementation of the polarization source term. 

(i) - "Linear source term". For the numerical results shown in Fig. 2 and in Fig. 7(c)-(d), the polarization source 
term is implemented as a linear relativistic inhomogeneity (RI), travelling at velocity v along the z direction, i.e., 

P{z,t) = eo{[no + 5n{t - z/v)]'^ - nl}E{z,t), (5) 

where the background term, no, is the refractive index at the input carrier frequency and t = t — z/v is the local 
time coordinate. The "artificial" RI (Sn) has a super-Gaussian longitudinal profile of amplitude Shq and thickness a: 
6n{T) — (5rioexp[— (r — <o)'"/c"']; being m a positive even number identifying the Gaussian order and to the initial 
temporal position. 

(ii) - "Nonlinear source term,". The numerical model is extended to inchide the nonlinear case presented in 
Fig. 7(a)-(b); that is the case in which the inhomogeneity is induced by an intense soliton through the nonlinear 
Kerr effect. The polarization source term is P{z, t) = eonon2\£{z, t)\'^E{z, t) where £ is the complex valued analytical 
signal and {l/2)n2\£(z,t)\'^ gives the dimensionless refractive index variation. 

The probe pulse scattering by the soliton-induced RI, is modeled by solving a second equation for the real field of the 
probe pulse, E2{z,t), coupled to the equation for the soliton (4): 

dzE2,u; - i[kz{oj) - u/v]E2,u; = o, / x ^2,u- (6) 

The coupling is indeed given by the "cross-phase modulation" term of the nonlinear polarization, i.e., P2{z,t) = 
2eonon2\£{z,t)\^E2{z,t), where £ is the complex valued analytical signal of the pump (soliton), as in the previous 
case. 

We intentionally neglect all the other nonlinear terms, e.g. third harmonic generation and Raman scattering, in 

order to highlight the generation of the RR and NRR modes. Moreover, we recall that, in order to capture the 
coupling between positive and negative frequencies, the model equation (and especially the source term), must be 
implemented by calculating the w-frequency component of the real polarization P{z,t). This is achieved for each 
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propagation step by means of back and forth Fourier transforms allowing for an evaluation of the real electric field 
E{z,t) from its complex spectral components Ef^{z). 
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ADDITIONAL MATERIAL 
A: Dispersion relations. 

As in all the scattering processes, the input and output scattered modes described in the main text, should obey 
precise momentum conservation laws. In the case of light scattered from a travelling inhomogeneity, the wave-vector 
k{ujji^) and frequency wrr of the scattered wave are linked to the input pulse by the condition: 

fc(wRR) = fc(t^iN) H , (^Ij 

V 

where k(uj) = u)n{uj)/c is the medium dispersion relation, v the relativistic inhomogeneity (RI) velocity, and win the 
input carrier frequency. The last term in Eq. (SI), (wrr — iOm)/v, is related to the momentum transferred to the 
incoming pulse from the RI [1]. This relation is rigorously derived in the following section. 

In order to understand the origin of the negative RR, it is useful to find its coordinates on the dispersion curve k{uj) 
by graphically solving Eq. (SI). Without any loss of generality, in Figure 5 we consider a simple dispersion relation 
such as that of diamond, which does not exhibit any resonances over a very broad bandwidth (from less than one 
Terahertz through to ultraviolet wavelengths). Equation (SI) is a straight line when represented in (w, k) coordinates, 
as in Fig. 5(a) (blue solid line). The frequency and wave- vector values of the propagating modes are thus determined 
by the overlap of this line with the actually allowed modes in the medium, i.e. by the intersection with the medium 
dispersion relation (red solid curve). Red dashed curves indicate the negative frequency branch (in the lab. frame) 
of the dispersion relation, while the grey shaded areas indicate the domains of the backward propagating modes 
(namely, the 2nd and 4th quadrants), which are not being excited by the forward propagating input mode and are 
not considered in this work. A first intersection is identified as the incoming mode IN. A second mode can be found, 
RR to which light may be resonantly scattered. Although the generation of resonant radiation is typically considered 
to be a property of self-scattering from intense nonlinear soliton pulses, we note that actually it is a more general 
process that only requires the existence of some kind of inhomogeneity in a linear medium, as discussed here. A third 
intersection is found on the negative frequency branch of the dispersion curve: this is the negative frequency resonant 
radiation mode, NRR*. Light may be resonantly transferred to this mode too, alongside the RR mode, as observed 
in recent experiments [2]. We stress that the existence of a solution at negative frequencies is not unphysical, being 
the electromagnetic field a reaZ-valued quantity, that naturally oscillates at positive and negative frequencies. Indeed, 
the mode that is actually measured in an experiment is obtained by adding the three modes described above with 
their complex conjugate counterparts. These are found from the intersections of the dispersion curve with the straight 
dashed line in Fig. 5(a), obtained from Eq. (SI) by substituting uj with —a; and using the fact that the dispersion 
relation is an odd function, i.e. /c(— w) = — fc(w). Physically, this describes the momentum conservation condition in 
which the complex conjugate of the IN mode, IN*, is being scattered towards RR* and NRR. 

In Fig. 5(b) we then consider the same dispersion relation in the reference frame that is comoving with the inho- 
mogeneity, i.e. in (w', fc') given by the Doppler transformations: lu' = 7(0; — vk) and fc' = 7(fc — ujv/c^), where 




FIG. 5: (a) Dispersion relation of diamond in the lab. frame (w, k) and (b) in the reference frame {uu' , k') moving at velocity v. 
Solid (dashed) red curves correspond to positive (negative) dispersion branch in the lab. frame. Straight blue lines indicate the 
momentum conservation condition, Eq. (SI) (solid blue), and its "phase-conjugate" counterpart (dashed blue). Grey shaded 
areas indicate the domains of backward propagating modes in the lab. frame, which are not excited by a forward propagating 
input pulse and are not considered in this work. 
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7 = 1/ — (w/c)2 and c is the speed of light in vacuum. In this reference frame the momentum conservation 
relation (SI) transforms into a condition on the comoving frequencies, w' = w'in, and therefore represented as a 
horizontal line in the graph. Specifically, we find w^j^ = +oo'i^ and = +w{j>^. The NRR mode thus satisfies 

'^NRR ^ ^'^IN- ^^'^ R-R- (NRR) mode is thus naturally described by a positive [negative) frequency in the comoving 
reference frame, reflecting the fact that RR (NRR) builds from frequencies with positive (negative) local phase ve- 
locities. We report here also the Sellmeier equation for the dispersion relation of diamond used in our simulations: 
n2(A) = 1 +4.3356 AV(A2-0.10602) +0.3306 AV(A2- 0.17502), where the wavelength A = 2-kc/(jO is expressed in jim. 



B: First Born approximation. 

In scattering theory the linear relation involving input and output modes may be modeled with the scattering 
matrix S, which is a unitary matrix connecting asymptotic states and which describes all the scattering channels of 
the process. Physically the square moduli jS'j.ip of the S-matrix elements represent the number of photons for the 
J— output mode, which is being scattered by the i— input mode. Following Kolesik et al. [3], we evaluate the elements 
of the S-matrix in the first Born approximation, giving an analytical formiilation for the spectral amplitude of the 
scattered wave. We also generalize the treatment in order to include the possibility of a third output scattered mode, 
i.e. the NRR mode that, as argued above may also be generated during the scattering process. 

The role of the "scattering potential", R{z,t), is played by the moving rclativistic inhomogeneity (RI), which we 
assume to be stationary along propagation, i.e. R{z,t) = R{t — z/v) = dn{t — z/v). 

In detail, let us consider the laser pulse evolution "unidirectional pulse propagation equation" (UPPE) in the labora- 
tory reference frame: 

where 

P{z,t) = eo {[nQ + 5n{t~ z/v)f- - nl) E{z,t) ~ 2eQno6n{t - z/v)E{z,t), (S3) 

so that 

P„ = -L f dt 2eono6n{t - z/v)E{z,t)e''^K (S4) 
2-K J 

Equation (S2) is the so-called "unidirectional pulse propagation equation" described in the Methods section and also 
numerically solved in this work. We want to solve this equation perturbatively by writing 

E^{z) = E^SHz) + £;W(^) + Ei'\z) + ..., (S5) 
where e2\z) is of order ((5n)'. We then get the equations: 

d^E^^^ -ih{oj)E(^^ =0, (S6) 

d,E(J)-zk,{oj)E(J^=i^^^^ I dt [<5n(t-^/^)e-*y'e— '*4°)(z)da;' 



(S7) 



Prom (S6) we get 



i?(o)(z)=e''^-(-)^i?o, (S8) 



where E^ is constant. Now, we can solve Eq. (S7) and, assuming that for z — >■ — cxd there is only the unperturbed 
solution, we find, 

Ej^\z) = i 2Sk\ujf '''^'^^^ f Sn{t-^/v)e'^* J e-'^''+"''^'^">^ E°, dtdcu'^ . (S9) 

We can rewrite this as 

E(;}\z) = Jdw'a{z,co,u;')Ej^}{z), (SIO) 
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where a is the S-matrix density: 

ZTTC kz\<jj) J J-oo 

By using the new set of variables, u = t — ^/v and w = t + ^/v, we can evaluate the integrals in the last expression, 
to obtain: 



(5 1 w - Win - V ■ [^^(w) - ^^(win)]! ^ |^ ^ ~ ^in + v ■ [kz{ui) - ^^(win)] | 



(S12) 



where S is the Dirac delta function, which is nonzero only for co — Win — v[kz{uj) — A;z(win)] = 0, and R is the Fourier 
transform of the scattering potential 5n{u), 

R{u) := / dw(5n(M)e*"". (S13) 



We can now substitute (S12) in Eq. (SIO) and by performing the integration (using the delta function) we finally 
obtain: 

E^^\z) = S{z,>^,u:,^)E(Z{z), (S14) 

where the S-matrix element is: 

S{z, w, win) = *^|^^e^^--™)-/"i?(w - win), (S15) 

and 

- — ^ = kz{uj) - kz{(jJm). (S16) 

V 

By substituting the delta function condition for w spanning over positive and negative values, we then retrieve 
Equations (1) and (2) of the main text. 

Summarizing, the scattering matrix elements that describe the scattering of a laser pulse from a RI, may be explicitly 
evaluated from Eq. (S15): the amplitude of the scattered wave is directly proportional to the Fourier transform of 
the scattering potential, i.e. of the RI. Moreover, Eq. (SI) emerges as natural condition that has to be satisfied in 
order for the scattering process to occur. This therefore provides a firm grounding for the fundamental momentum 
conservation relation used in this work, Eq. (SI). 

C: Generalized Manley-Rowe relation. 

For the scattering process under study, starting from an input seed pulse at frequency win, the photon number of 
the output waves is linked to the photon number of the input mode, |INp, by the following Manley-Rowe relation, 
generalized to the case of a moving scatterer [4, 5] : 

sign(7RR)|RR|' + sign(7NRR)|NRR|2 = sign(7iN), (S17) 

where |RRp and |NRRp arc the photon numbers of the RR and NRR modes normalized to the input photon number, 
|INp, and 7i = (w| — vv^)i, where is the phase velocity of the i-mode. We note that sign(7i) = sign(wi — vki) = 
sign(w^); while the real- valued RR mode will have positive frequency in the comoving frame, i.e. w^j^ > 0, the 
real-valued NRR mode will always have a negative comoving frequency, w^rr < 0. Therefore, Eq. (S17) assumes the 
form of Eq. (3), that is [RRp - [NRRp = 1. 
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FIG. 6: Propagation cross-correlation maps for (a) linear scattering in diamond and (b) nonlinear scattering in fused silica. 
The two maps refer to the numerics presented in the main paper, Figures 2 and 3(a)-(b), respectively. Arrows indicate the 
position of the IN, RR and NRR modes. 



D: Propagation Cross-Correlation. 

We evaluate the propagation cross-correlation map for the numerics presented in the main text, by adopting the 
approach of Bejot et al. in [6]. This correlation is positive if frequencies are created simultaneously, i.e. at the same 
propagation distance z, and negative if one frequency is used to create the other. For each pair of frequencies the 
propagation cross-correlation is defined as: 

C(.„..)= , Cov{/K),/(..)} 

v/Cov {I{uj,),I{oj,)} Gov {I{oj^),I{u,)} 

where, 

Cov{/(c.i),/(c.2)} = {[liu^uz) - {I{oJi,z)}] [I{oj2,z) - {I{i02,zm, (S19) 

where (•) indicates the expectation value i.e. the average evaluated over the propagation distance z Figure 6(a) shows 
the propagation cross-correlation map for the numerical results reported in Fig. 2 of the main paper, which refer to a 
linear scattering in diamond. The two resonant modes, RR and NRR, are perfectly correlated (correlation equal to 
-f 1), implying that they are born together, whereas they are both anti-correlated with the input mode. This confirms 
that indeed the RR and NRR modes are the result of the same scattering event seeded by the IN mode. Finally, 
Figure 6(b) shows the propagation cross-correlation map for the numerics of Fig. 3(a)-(b) of the main paper, referring 
to a nonlinear scattering in fused silica. Here too, the map shows obvious signs of correlation between the resonant 
modes. 



E: Comparison between a soliton-induced and a linear, "toy- model", RI. 

Figure 7 is an extension of Figure 3 of the main paper and shows a comparison between the dynamics of a soliton- 
induced RI, with self-scattering towards the resonant modes (a)-(b), and an artificial, "toy- model" RI co-propagating 
with a weak probe pulse in a purely linear medium (c)-(d). Details of the simulations are indicated in the main text. 

For both cases the frequencies of the RR and NRR modes are in perfect agreement with the theoretical predictions, 
as shown on the conioving dispersion curve (a;,^') in Fig. 7(f). Finally, Figure 7(e) [as Fig. 6(b)] shows the propagation 
cross-correlation for the nonlinear case, as described in detail in Section D and in [6]. The two resonant modes are 
correlated (correlation equal to +1) implying that they are born together, whereas they are both anti-correlated 
(correlation equal or close to —1) with the input mode. This provides further confirmation that indeed the RR and 
NRR modes are the result of the same scattering process seeded by the IN mode. 
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FIG. 7: RI scattering dynamics in fused silica. Comparison between a soliton-induced (via the nonlinear Kerr effect) RI 
(a)-(b) and the "toy-model" RI with the same shape as the soliton-induced inhomogeneity, in a purely linear medium (c)-(d). 
Envelope (left hand side), and spectral (right hand side) evolution in propagation. Propagation cross-correlation map for the 
nonlinear case (e). Theoretical positions of input (IN) and resonant (RR, NRR) modes on the comoving dispersion curve 
oj' = uj' (uj), indicated with filled circles (f). 
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